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ABSTRACT 


Equations  are  derived  for  predicting  the  decay  day  and  decay 
revolution  of  a  satellite.  Only  two  parameters,  obtained  from  the 
quadratic  equation  that  predicts  equatorial  crossings,  are  UBed  . 
Examples  of  predictions  are  given. 
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1.  INTRODUCTION 


<4 


* 


Equations  are  derived  for  predicting  the  decay  day  and  decay 
revolution  of  a  satellite.  Only  two  parameters,  obtained  from  a 
quadratic  equation  which  predicts  equatorial  crossings,  are  used. 
It  appears  that  a  predicted  decay  day  might  be  useful  in  predicting 
equatorial  crossings.  Some  examples  of  predictions  are  given, 

2.  DERIVATION  OF  EQUATIONS 


Let 


.PN  =  F{TN)  (1) 

where  F  (Tpj)  is  a  continuous  function  of  time  and  pN  is  the  time 
derivative  of  the  nodal  period  at  revolution  number  N.  Integration  of 
(l)1,  by  the  first  law  of  the  mean  for  integrals,  yields 


where 


•  TN  =  PN  '  PN  =  PN  -  PN  » 


i  d  i 
F(TN  ) 
j 


i  d 


(2) 


E  =  exact  number  of  days  to  decay  from  revolution 
Ni 

number  N^, 

Nj  =  revolution  number, 

N(j  =  revolution  number  at  decay, 

=  time  at  decay  (in  days), 

1  Figures  in  parentheses  refer  to  equations. 


I 


T  =  time  at  revolution  N  (in  days), 


N, 


PM  =  nodal  period  at  TM  (in  days). 


N, - -  ** - -  *  N. 

1  i 


P  =  period  at  TN  (in  days), 
wd  d 


Pfl  =  mean  value  of  during  the  interval  E  ,  and 
j  N  Nj 


Set 


LN,  *  TN  *  V 
d  j  i 


PN.  (PN.  ■  PN  ) 
d  -  i  i  d 

PN,  “ - 


'j 


°-6<1+aNi)  *PN.  “  lffo 


<3) 


where  P.^  is  the  rate  of  change  of  the  period  at  TM  ,  and 


N. 


aNt  =  aNi  ^PNj  *  PNj*  PNd>-  <4) 

/  88 

The  constants  in  (3)  tend  to  minimize  the  OC^  when  P^.  iB  greater  than  days 

Substitution  of  (3)  into  (2)  gives 


88 


i^N.  *  0.6  (1  +  ajyj _)  (Pjsjj,  -  1440  ) 


(5) 


-  P 


Nj 


For  predicting  the  time  of  equatorial  crossing,  The  National 
Space  Surveillance  Control  Center  (NSSCC),  (L.  G.  Hanscom  Field, 
Bedford,  Massachusetts)  uses  the  equation 

tn  +  n  =tn  +  Npn  +  n2  cn  +n3  dn  .  (6) 

i  i  i  i  i 


2 


where  is  generally  omitted  until  the  last  hundred  or  so  revolutions 
o£  a  satellite. 

Define  in  (6)  by 


1 


Define  Aj^  by  either 


88 


0,6  *PN-  “  1440  > 


kN, 


i 


-Pw 


or 


An  . 


"N 


Substitution  of  (8)  ir.to  (5)  gives 


(8) 


(9) 


ENi  “  +  “Nj^  ANa  •  <10> 

No  assumed  value  of  decay  is  contained  in  (10). 

In  attempting  to  obtain  an  empirical  representation  for  , 

through  an  analysis  of  data  from  15  decayed  satellites,  the  predictions 

A  88 

indicated  that  |oj^  |  <  if  the  predictions  were  made  when  Pj^  >  -^q 

days.  When  this  condition  is  fulfilled,  the  accuracy  of  the  predictions, 
given  by  (9)  indicates  that  (9)  is  a  good  approximation  to  the  exact  solu¬ 
tion  (2)  for  those  15  satellites. 

Assume  for  a  given  satellite  the  P^  and  Cj^  in  (11)  are  such  that 
predicted  equatorial  crossings  are  in  good  agreement  with  observations. 
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*N  +  H,  “  TNt  +  N  PN:  ♦  N‘  CN; 


Substitution  of  these  values  of  and  C  into  (9)  gives  the 

Ni  Ni 

days  to  decay  from  T.r  .  The  decay  day  T._  is  then 
Ni  Nd 


TNd  =  TNi  +  ANt  ‘ 


The  Tj^  in  (12)  can  be  used  to  determine  Cj^  at  some  later  epoch 
Tj>j  ,  From  the  NSSCC  bulletins  we  obtain  TN  and  Pj^  .  Since  Tj^ 


is  now  known,  we  have 


Awk  =  TNd  ‘  TNk 


r„  -  “S'X'B'1  ’  (14) 

TNk  -  ^Nd 

Appendix  II  gives  an  application  of  this  technique  for  determining  the 
C«t  to  be  used  on  later  bulletins. 


3.  THE  DECAY  REVOLUTION 

In  order  to  predict  the  approximate  decay  revolution  number  N^, 
use  (12)  to  obtain  T,.  .  Then  let  N  be  the  revolutions  from  Nj  to 

a  ^ 

Nd  and  solve  (15)  for  Nq  . 

TNd  =  +  Nq  =  TN.  +  Nq  PNi  +  Nq*  CN.  .  (15) 
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(16) 


« 


t 


* 


Nd  =  Nq  +  Ni  • 

For  a  better  approximation  of  Nd,  especially  when  predicting 
far  in  advance  of  decay,  the  above  method  should  be  modified  in  the 
following  manner.  We  first  derive  an  expression  for  in  terms 


°f  Pj»j^  and  Cjyj  .  This  derivation  requires  that  the  Pj^  *n  ^7)  be 


88 


assigned  the  value  Y^q"-  days.  When  this  is  done,  the  approximate 
solution  (9)  is  exactly  the  same  as  that  obtained  if  (1)  iB  replaced  by 

Pn  =J  <PN“  PN  >  ‘  2  3  *  <17> 


where  J  is  a  constant.  From  (17)  we  obtain 

P  _  ip  _§«_  v"1  ip  ? 

PNi  “  _ ^  'PN^  ‘  1440  >  lpN^ 


(18) 


Define  D  in  (6)  by 
Ni 


PN- 


(19) 


Substitution  of  (7)  and  (13)  into  (19)  gives 
—  4  2  in  88  »  *  1 


D 


Nj=  -f  cN.«  (PNl-T^o> 


(20) 


We  substitute  D„.  into  (21)  and  solve  for  N 
Ni  q 


XNd  11  TNq+  N.  -  TNj  +  Nq  PNs  +Nq  CNi 


+  V  DN.  • 

J 


(21) 
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This  when  substituted  into  (16)  will  in  general  give  a  better  approxi¬ 
mation  of  than  will  the  Nq  obtained  from  (15).  To  illustrate  the 
last  statement  v/e  use  1959  Epsilon- 2  as  an  example.  The  first  two 

p 

NSSCC  bulletins  for  this  satellite  gives  for  epoch  revolution  number 

2, 

T2  =  46.08599 

P2  =  7.241102  x  10’2  (22) 

C2  =  -  6.795  x  10'7  . 

Substitution  of  (22)  into  (9)  gives 
A2  =  361.25278  days. 

Equation  (IS)  becomes 

361.25278  =  Nq  P2  +  Nq*  C2  (23)  , 

From  (23)  we  obtain  Nq  =  5247,  and  from  (22)  and  (20)  we  obtain 
D2  ®  -1.816  x  10  _11. 

Substitution  of  Dg  into  (21)  gives  Nq  a  5289.  Hence  the  first 

method  gives  N^  =  5249  and  the  second  method  gives  Nd  -  5291  . 

Since  the  decay  revolution  was  approximately  5313.5,  it  is  quite 

apparent  that  the  introduction  of  D2  gives  a  better  prediction  of  N^. 

It  should  be  pointed  out  that  (20)  may  also  serve  a  useful  purpose 

in  determining  the  values  of  D  to  use  in  (6)  when  predicting  equatorial 

Ni  * 

crossings  during  the  last  hundred  or  so  revolutions  of  a  satellite. 
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4.  REMARKS 


Some  of  the  variations  in  the  predicted  decay  day,  (obtained  by 
using  ,  PN  ,  and  C«r  from  NSSCC  bulletins),  are  due  to  the  fact 

IMi  j  lNi 

that  bulletins  must  be  issued  even  though  insufficient  or  no  observations 
are  available  on  which  to  determine  good  values  of  ,  Pj^,  and  Cj^. 

Since  the  predicted  decay  day  depends  only  on  these  values,  the  best 
available  values  should  be  used.  Such  values  are  obtained  by  fitting 
(11)  to  the  observations.  An  example  of  this  is  the  prediction  of  decay 
for  1958  Epsilon  in  Appendix  I. 
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APPENDIX  I  -  EXAMPLES  OF  DECAY  PREDICTIONS 
1958  EPSILON  (EXPLORER  IV) 

Equatorial  crossing  times  for  the  first  937  revolutions  of  this 
satellite  are  contained  in  "Interim  Definitive  Orbit  for  the  Satellite 
1958  Epsilon  IV  "  NASA  TN  D-410.  From  this  reference  we  obtain 
for  epoch  revolution  100,  Tj^q  =  3. 267095  August  1958.  The  equation 

TN  +  100  =  T100  +  N  P100  +  N*  C100 
predicts  T200  z«ro  error  and  T3QQ  with  an  error  of  14  seconds, 
when  P100  a  7.63746  x  10'2  and  Cjqq  =  -  7.84  x  IQ-7  .  Substitution 
of  these  values  into  (9)  gives  A10q  ■  446.1  (dayB  to  decay  from  Tjqq)> 
The  decay  day  was  approximately  23,2  October  1959,  so  the  error  in 
A100  iB  le8S  than  0.4  days. 

The  following  predictions  were  all  obtained  by  using  the  TM  , 

Pjq^  •  and  CN^  given  on  NSSCC  bulletins  issued  at  epoch  revolution 
number  N^  . 

1958  DELTA  (Sputnik  III) 


Ni 

an. 

1 

XNd 

7400 

171.68 

463.17 

7600 

147.66 

452.52 

8100 

131.07 

469. 12 

8250 

115.09 

463.03 

8 


\ 


t 


8700 

84,15 

461,54 

8900 

71,88 

462.24 

9000 

66.19 

463.01 

9200 

53,57 

463.26 

The  average  value  of  is  462,2.  The  actual  decay  of  the  satellite 
occurred  approximately  at  462,3. 

1959  ZETA  (Discoverer  VI) 


Ni 

ANj 

TNd 

30 

51.81 

285.57 

80 

59.11 

296. 18 

130 

57.05 

297.41 

190 

55.31 

299.61 

340 

46.65 

300.75 

500 

27.54 

292.00 

610 

17.  54 

289.06 

700 

13.84 

291.09 

720 

14.01 

291.39 

The  average  value  of  T  is  293.67.  The  last  NSSCC  bulletin  esti- 

Nd 

mated  the  decay  occurred  about  293.75. 
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1959  KAPPA  (Discoverer  VII) 


N. 

l 

an. 

tn. 

i 

d 

lb 

22.73 

335.54 

30 

11.79 

325.59 

50 

15.59 

330.69 

110 

12.19 

331.84 

200 

5.84 

330.61 

The  average  value  of  is  330.85.  The  last  NSSCC  bulletin  estimated 
that  decay  occurred  shortly  after  330.78. 


1959  EPSILON- 2  (Capaule  of  Discoverer  V) 


In  Appendix  II  we  find  that  the  first  two  bulletins  predict  A2  =  361.25 
and  TNd  =  407.34,  Since  these  bulletins  gave  accurate  predictions  on 
equatorial  crossings,  assume  that  the  decay  day  is  (407  +  30).  The 


average  value  of  Tj^  obtained  from  the  first  12  bulletins,  which  predict 
within  the  above  limits,  is  408,33.  All  twelve  predictions  were  made 


over  225  days  before  decay.  The  actual  decay  was  T-,  =  408.77. 

wd 
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APPENDIX  II  -  EXAMPLE  OF  PREDICTED  DECAY  DAY  USED  TO 
PREDICT  EQUATORIAL  CROSSINGS 

9  In  this  appendix  we  illustrate  how  the  predicted  decay  day  can 

be  used  to  obtain  the  CNi  terms  for  later  bulletins.  We  use  1959 
Epsilon-2  (Capsule  of  Discoverer  V)  as  an  example.  The  first  revolu¬ 
tion  of  this  satellite  was  arbitrarily  chosen  as  the  first  revolution  on 
15  February  I960.  The  first  and  Becond  bulletins  issued  by  NSSCC 
gave  good  predictions  over  an  interval  of  173  revolutions.  Both 
bulletins  used  T2  =  46.08599,  P2  =  7.241102  x  10"2.  and  C2  = 

-  6,795  x  110'7.  From  this  data  we  compute  A2  =  361.25298  and 
Tp^j  =  407.33897.  From  bulletin  #3,  we  obtain  Tj^q  =  56,06614  and 
*  P140  =»  7-22269  x  10‘2.  Then 

„  0.3  Pi40  (pi40  -  jJJj-) 

u140  - - - - 

T140  '  XN, 

Q 

=  -  6.06  x  10*7  . 

We  use  (11)  which  becomes 

T200  =  T140  +  60  P140  +  3600  C140, 
to  predict  T2qq’  °n  bulletin  #5,  T20Q  =  71,19978.  The  error  in  the 

prediction  is  21  seconds, 

s 

On  bulletin  #33,  Tjqqq  =  256.91272  and  Pjqqq  =  x  10  2  . 
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We  use  the  above  and  repeat  the  procedure  to  predict  the  epoch 

on  bullet1  n  #35,  T310Q  =  263.73626.  The  error  in  the  prediction  is  2 
seconds. 

On  bulletin  #48,  T  . ,  =  330.76048  and  P  a  6. 65805  x  10'2. 

4100  4100 

We  repeat  the  above  procedure  to  predict  the  epoch  T4200  “  337.32994 
on  bulletin  #S0,  The  error  in  the  prediction  is  24  seconds. 

An  application  of  the  above  to  four  satellites  Btill  in  orbit  indicates 
that  this  technique  may  have  some  value. 
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APPENDIX  IM  -  DECAY  PREDICTION  PROGRAM  WITH  EXAMPLES 
In  this  appendix  we  discuss  a  program  for  obtaining  three  pre¬ 
dictions  of  Tn  from  one  epoch  TM  .  BN  is  the  TN  predicted  by 
iNd  Ni  ”i  iNd 

the  Tm  >  PM  ,  Cm  given  on  the  bulletin.  B..  is  the  Tv,  predicted 
Ni  Ni  Ni  Nt  Nd 

if  CM  is  changed  to  CN  in  order  that  (11)  will  predict  TM  with 
wi  i  Ni  +  1 

zero  error.  is  the  TM  predicted  if  Cm  is  changed  to  Cm  in  orde 

Ni  wd  i  i 

that  (11)  will  predict  T  with  zero  error. 

Ni  +  2 

Let  P  ,  T  ,  C  ,  (j  =  0,  1,  2),  be  the  elements  on 

Ni+j  Ni+j  Ni+j 

three  successive  bulletins. 


Compute 


=  -3Pm  (pNt~ nsr*--0*,  ^n,)'1  • 


BNi  =  +  AN.  » 


_  tN.  .  -  tN.  -  (Ni  +  i  -  Nt)  PN  . 

^  _  1  +  1  1  1 
^  “  “ 


Mi  +  1  •  "iJ 


AN'  =  ’  °N- 

1  1 


bn.  =  tn.  +  an.  * 
1  1  1 
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c. 


\  t  ,  -  Tn,  -  (N1  *  »  •  ^ 


N, 


(Ni  +  2  - 

ANj  =  '  GNi  (CNt)  1  *  and 


BNt  "  T»i  +  Aui 


The  equations  were  programmed  for  the  IBM  709  computer.  The 
printout  gives  the  above  eight  quantities  and  also  the  orbital  elements 
contained  on  the  bulletins.  The  program  starts  with  1=1,  where  Nj 
is  the  epoch  revolution  on  the  first  bulletin,  and  then  repeats  the  pro¬ 
cedure  for  all  values  of  i. 

It  was  believed  that  the  program  might  provide  B  and  B  values  which 
gave  a  more  accurate  predicted  decay  day  than  did  the  B  values. 

Additional  calculations  are  necessary  before  the  program  can  be  evaluated. 
Some  predictions  calculated  by  this  program  are  given  below. 

1959  LAMBDA  (Discoverer  VIII) 


A780  =  39,8 

B780=54-3 

X780  =  53.  8 

"®780  =  68-3 

A780  =  51 1 7 

B780  “  86,2 

Decay  day  was  approximately  68.0. 


> 

.1 


) 
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1959  EPSILON- 1  (Discoverer  V) 

A9  =  29.1 

B9  »  255.4 

"X9  b  45. 5 

•ET0  «  271.8 

A9  »  44.7 

B9  b  271.1 

Decay  day  was  between  271  and  272, 

1958  EPSILON  (Exolorer  IV1 

A4170  =  150,0 

B4170b  297-4 

"*4170  =  141,1 

^4170  B  288*5 

A4170  *  150,6 

b4170  =  2^8,0 

Decay  day  was  approximately  296.2. 
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